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Abstract —When dedicated positioning systems, such as GPS, 
are unavailable, a mobile device has no choice but to fall 
back on its cellular network for localization. Due to random 
variations in the channel conditions to its surrounding base 
stations (BS), the mobile device is likely to face a mix of both 
favorable and unfavorable geometries for localization. Analytical 
studies of localization performance (e.g., using the Cramer-Rao 
lower bound) usually require that one fix the BS geometry, and 
favorable geometries have always been the preferred choice in 
the literature. However, not only are the resulting analytical 
results constrained to the selected geometry, this practice is likely 
to lead to overly-optimistic expectations of typical localization 
performance. Ideally, localization performance should be studied 
across all possible geometric setups, thereby also removing any 
selection bias. This, however, is known to be hard and has 
been carried out only in simulation. In this paper, we develop 
a new tractable approach where we endow the BS locations 
with a distribution by modeling them as a Poisson point process 
(PPP), and use tools from stochastic geometry to obtain easy- 
to-use expressions for key performance metrics. In particular, 
we focus on the probability of detecting some minimum number 
of BSs, which is shown to be closely coupled with a network 
operator’s ability to obtain satisfactory localization performance 
(e.g., meet FCC E911 requirements). This metric is indifferent to 
the localization technique (e.g., TOA, TDOA, AOA, or hybrids 
thereof), though different techniques will presumably lead to 
different BS hearability requirements. In order to mitigate 
excessive interference due to the presence of dominant interferers 
in the form of other BSs, we incorporate both BS coordination 
and frequency reuse in the proposed framework and quantify 
the resulting performance gains analytically. 

Index Terms —Cellular localization, E911, hearability, stochas¬ 
tic geometry, point process theory, base station coordination, 
frequency reuse. 

I. Introduction 

G EOLOCATION (also called positioning, localization, and 
position location) is deeply ingrained in our daily lives 
and has been studied by the scientific community for many 
years @-© ■ The driving force behind much of the research is 
a mandate by the Federal Communications Commission (FCC) 
requiring cellular network operators to locate those calling 911 
to within certain accuracy requirements 0,0. Until recently, 
these requirements included only outdoor location accuracies. 
Accordingly, the predominant way cellular network operators 
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have met the requirements of the mandate is by relying on the 
Global Positioning System (GPS). In January of 2015, however, 
the FCC expanded its mandate to include a phase-in of indoor 
positioning requirements, citing that the bulk of emergency 
calls now originate indoors |9|]. 

While accurate outdoor positioning using GPS is reliably 
available under clear sky conditions, many years of positioning 
study have not yet resulted in equally reliable positioning 
methods in urban canyons and indoor scenarios. Consider 
the classical example of locating emergency personnel (e.g., 
firefighters) indoors, where our current inability to provide 
accurate indoor positioning has recently been described as a 
dilemma “where people [are] literally dying within a hundred 
feet of safety’ ’G3- Despite this, economic limitations dictate 
that global navigation satellite systems (GNSS), such as 
GPS, are likely to be the only widespread dedicated location 
systems in the foreseeable future. This necessitates a fallback 
to terrestrial cellular networks for geolocation in situations 
where these prevalent location technologies are unavailable. 
Presently, however, no analytical approaches exist to study the 
fundamentals of localization performance in these networks. 

It is the objective of this paper to introduce a new 
tractable model for studying terrestrial geolocation using 
cellular networks. The model uses concepts from point process 
theory E) and stochastic geometry [ [12] to lend tractability 
to the study of network localization, which is demonstrated 
through accurate and easy-to-use expressions characterizing the 
number of base stations (BSs) that are able to participate in a 
localization procedure, a key factor influencing the localization 
performance. 

A. Prior art and motivation 

Regardless of the technique, geolocation performance fun¬ 
damentally depends upon three things: (i) the number of 
participating BSs, (ii) the geometry or locations of these BSs 
relative to the device being localized, and (Hi) the accuracy of 
the positioning observations. When these factors are determin¬ 
istic, localization performance is typically studied analytically 
using the Cramer-Rao lower bound (CRLB) (e.g. G3-(B])- 
Furthermore, since localization systems are being discussed, 
the topologies considered are often ones favorable to location 
estimation (such as by using a carefully controlled set of 
scenarios to avoid poor geometries, e.g. E))- This approach, 
however, is not appropriate for studying cellular localization 
performance because none of the aforementioned factors may 
reasonably be considered deterministic. Instead, the different 
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possible locations of the mobile device within a network, the 
topology of the BSs relative to the mobile device, and random 
variations in the channel conditions will (i) affect the number 
of BSs capable of participating in the localization procedure, 
(ii) result in those BSs providing a mix of both good and bad 
topologies for localization, and (iii) define the strengths of the 
received BS signals which directly impact the accuracies of 
the positioning observations. Clearly, only a very large set of 
deterministic scenarios could provide insight into localization 
performance over all possible variations of the above. Even so, 
crafting realistic scenarios requires the assumption of specific 
system design parameters and propagation conditions, meaning 
that any insights from the results would be specific to the 
chosen assumptions, rather than being general. This makes a 
classic metric such as the CRLB not suitable for a general 
analysis. 

In this work, our focus is on providing a truly general analysis 
of one metric which we will see is indicative of a network 
operator’s ability to meet its mandated localization accuracy 
require me nts-the number of successfully detectable BSs. For 
communication, it is ideal for a mobile device to receive a 
strong signal from its serving BS and weak signals from all its 
neighbors. For localization, on the other hand, it is increasingly 
beneficial for a mobile device to receive usable signals from 
more and more neighbors. Although communication is the 
primary concern of network operators, cellular system designers 
have no choice but to also cater to localization demands. As a 
result, the need to satisfy these conflicting goals has a name 
among designers-the hearability problem , due to the fact that 
increasing hearability from neighboring cells for the purposes 
of positioning is contrary to the principles of cellular system 
design In order to provide a general analysis of cellular 
hearability, it is appropriate to characterize this metric across 
all possible network topologies and channel conditions for the 
aforementioned reasons. 

One field of study in which random network topologies are 
considered when discussing localization is that of wireless or 
ad-hoc sensor networks. However, these studies differ from 
the cellular scenario we consider in that the sensor network 
literature almost certainly ignores interference and propagation 
effects, instead opting to define coverage using some fixed 
detection range 0, 0, (18), though in rare cases the 
received signal strength may be used to define a circular 
coverage region when shadowing is also considered 03 - 
In order to apply to cellular networks, our analysis differs 
by not only including interference, but also by employing a 
model enriched with additional cellular and localization-specific 
parameters. 

Lastly, when cellular networks are specifically considered, 
they are typically studied using the popular hexagonal grid 
model, making analytical results difficult to come by. Instead, 
network designers resort to complex system-level simulations 
using a common set of agreed-upon evaluation parameters 


as is clear from the 3GPP standardization process |2l |. This 
can also be seen in the limited amount of literature on 
hearability |22|-[24 1. Though simulation-based, the work in 
[ 22] and [|23| is most similar to ours in that the focus is on 
studying cellular hearability (in GSM and WCDMA (22) and 


LTE (23)). Interference, coordination, and fractional load are 
considered, but in the end, one is left with only simulation 
results, from which it is not possible to gather general insights 
into how localization performance is impacted by system design 
parameters and propagation effects. This motivates the need for 
a tractable analytical model that can provide preliminary design 
insights and will either circumvent the need for simulations 
completely or limit the ranges of the simulation parameters. 
In [[25J and [26), such tractable, yet accurate, analytical models 
were developed to study the coverage and rate of single-tier and 
multi-tier cellular communication systems, respectively. It was 
shown that the spatial layout of the BSs in cellular networks 
can be tractably and realistically modeled using Poisson point 
processes (PPPs), especially as cellular networks continue to 
deviate from centrally-planned macro-cell networks to networks 
which include an increasing number of more arbitrary small¬ 
cell deployments such as picocells and femtocells. This opens 
the door to the use of powerful tools from stochastic geometry 
to derive closed form expressions for key performance metrics. 
Motivated by these advances, we develop a similar approach 
to lend tractability to the analysis of localization systems. 


B. Contributions 

The main contributions of this paper are as follows. 

A general model for studying localization in cellular networks: 
In Section [TT] we build on the foundations of (25] and [ |26 ] and 
present a tractable model specifically for studying localization 
in cellular networks. The model includes average network load, 
the ability for BSs participating in the localization procedure 
to coordinate transmissions (e.g., through joint scheduling in 
LTE (16|), and is easily adapted to frequency reuse when 
frequency bands are assigned randomly. From |25|, it is easy 
to see that self-interference from the network often leads to 
poor coverage probabilities even when the reception of only a 
single signal is required. Since localization procedures require 
the successful reception of multiple signals and are thus even 
more demanding, localization systems typically necessitate 
the integration of signals over time in order to improve 
detectability. Our model assumes that this integration provides 
some fixed processing gain, while fast fading is excluded under 
the assumption that it is averaged out at the receiver. This 
exclusion of small-scale fading effects in our model agrees 
with the state-of-the-art 3GPP simulation models (20). 
Definition and characterization of accuracy-related hearability 
metric: We define a metric in order to study the number of 
BSs whose signals arrive with sufficient quality to successfully 
participate in a localization procedure. After presenting a clear 
relationship between the metric and positioning accuracy, upper 
bounds and approximations for the distribution of this metric are 
derived using the proposed model, resulting in expressions that 
are closed-form in some special cases of interest, while easy- 
to-compute integral expressions for the others. A dominant- 
interferer type analysis is employed (27), (28), whereby the 
strongest interference source is always treated exactly, leading 
to remarkably accurate approximations, nearly indistinguishable 
from truth. 

Design insights: Lastly, we present results and draw out some 
interesting observations useful in the design of localization 
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TABLE I 

Summary of Key Notation 


Notation Description 


|| ■ ||; a 

A 

o 

Xi (Xi) 

Ri ( Ri ) 
7 

P 

L 

P 

q 

n 

b(0,r) 
|.4|; A c 
BCA 
A\B 
1(.) 


£2 -norm; path loss exponent (a > 2) 

PPP of BS locations; Density of $ 

Origin (location of the typical user) 

Location of the ith closest (active) BS in $ 

Distance of x, (xi) from the origin 
Processing gain 

Post-processing SINR required for a successful wireless link 

Number of participating BSs 

Activity factor among L participating BSs 

Average network load for non-participating BSs 

Number of participating BSs interfering with the I A BS 

A ball centered at 6 with radius r 

The Lebesgue measure of region A', the complement of A 
Region B is contained in A 

Region A excluding its overlap with B\ A\B = A n B c 
Indicator function; l(-4) = 1 if A is true, 0 otherwise 


systems. We make evident the fact that localization performance 
is limited by the interference from the BSs participating in 
the localization procedure. It is shown that mitigating this 
interference through BS cooperation helps up to a point, beyond 
which the true gains in localization performance are realized 
through frequency reuse. The need to employ frequency reuse 
for accurate positioning has been observed previously, but only 
by first going through the process of running many complex 
system-level simulations. In LTE, for example, a reuse of six 
was deemed necessary for downlink TDOA positioning, a 
verdict reached after many simulations IZD- 

ii. System Model 

We now formally describe the proposed model. The key 
notation presented in this section is summarized in Table [I] 


A. Spatial base station layout 

The locations of the BSs are modeled using a homogeneous 
PPP $ £ R 2 with density A 1121. Due to the stationarity of a 
homogeneous PPP, the device to be localized is assumed to be 
located at the origin o. If the interference is treated as noise 
at the receiver, the most appropriate metric that captures link 
quality is the signal-to-interference-plus-noise ratio (SINR). 
For the link from some BS x £ $ to the origin, the SINR can 
be expressed as: 


SINR X 


-P£r||z||- a 

Eyet PSy\\y\\~ a + cr 2 ’ 


(1) 


y^x 

where P is the transmit power, S z denotes the independent 
shadowing affecting the signal from BS z to the origin, a > 2 
is the pathloss exponent, and er 2 is the noise variance. Note 
that 0 represents the SINR prior to any processing gain, yet 
as will be evident in the sequel, positioning systems typically 
have to work at lower target SINRs, thereby necessitating the 
need for some form of processing gain. In general, the post¬ 
processing SINR will include some multiplicative factor 7 
representing the processing gain, which depends upon system 
parameters (e.g. integration time) and is assumed to average 
out the effect of small scale fading. This is the reason why the 


SINR expression in (JTJ does not contain a fast fading term, 
which is consistent with current models for evaluating cellular 
positioning performance |20[ . Those conversant with stochastic 
geometry-based analyses of wireless networks will recognize 
that the absence of fast fading on the serving link, specifically 
one from the exponential family of distributions (e.g., Rayleigh 
fading), adds to the technical challenge of a localization system 
analysis. Lastly, although we consider interference as noise, 
the gain from some multiuser detection technique, such as 
interference cancellation, can be abstracted into the proposed 
model using the processing gain parameter 7 J29). 


B. Selection of participating base stations 

As is common in communication system analyses, serving 
BSs are selected according to the strongest BS association 
policy, measured using average signal strength, which typically 
includes long time-scale effects such as shadowing and pathloss. 
Now, consider that we desire to use a total of L BSs for 
positioning. Thus, we assume that the L BSs which provide 
the highest average received power make up the set of 
participating BSs. Their successful participation, however, is 
not guaranteed as there is some post-processing SINR threshold 
P (or equivalently, pre-processing SINR threshold /?/ 7 ) above 
which the signals from the participating BSs must arrive in 
order for them to successfully contribute to the localization 
procedure. In the absence of shadowing, the set of potential 
BSs simply corresponds to the set of the L nearest BSs. 
When shadowing is considered and BSs are selected according 
to average signal strength, the effect of shadowing may be 
absorbed as a perturbation in the locations of the BSs provided 


that the fractional moment E 


5 ; 


.2/0 


< 00 1 30 


when this condition is fulfilled and without loss 


of 
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. Thus, 
generality. 
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2/c 


. In 


we define a new equivalent PPP with density A E 
doing so, we ensure that the strongest BS association policy 
in the original PPP is equivalent to the nearest BS association 
policy in the transformed PPP. For notational simplicity, we 
will continue to represent the transformed PPP by $ with 
density A, under the assumption that if shadowing is present, it 
is already reflected in the density of $. Note that the condition 
on the fractional moment is fairly mild and will almost always 
be satisfied by the distributions of interest, including the most 
common assumption of log-normal shadowing with finite mean 
and standard deviation j3TJ. As a final note, it has been shown 
that for SINR, shadowing causes even more regular network 
models, such as the common hexagonal lattice, to behave like 
a PPP model |32]-|34). This further validates the use of a PPP 
to model BS locations. 


C. Base station coordination and network load 

BS coordination and network load are modeled through two 
activity factors, p and q, respectively. During a localization 
procedure, the L participating BSs coordinate with each other 
and attempt to blank their own transmissions while the others 
are active, but are unable to do so with probability p due 
to network traffic demands. The remaining BSs are assumed 
active with probability q, which is simply the average network 
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Fig. 1. The proposed model. The mobile device is denoted by the filled 
square and actively transmitting BSs are denoted by solid circles. 


load (35| . For simplicity of exposition, let us order the BSs 
in the now shadowing-transformed $ in terms of increasing 
distance from the origin such that the location of the k th farthest 
BS from the origin is denoted by Xk £ ‘l 1 - We now enrich the 
previous SINR expression in ([!]) for the signals arriving from 
the participating BSs (i.e., Xk for k £ {1,..., L }) by including 
BS coordination and network load as follows: 

SINR k(L) = - - P|l ^ fc11 " -, (2) 

J2a i P\\x i \\- a + bjPWxjl]- 01 + a 2 
J ~ L+1 

where ai and bj are independent indicator random variables 
(fixed throughout the localization procedure) modeling the 
thinning of the interference field by taking on value 1 with 
probabilities p and q, respectively. Note that SINR;, is now 
a function of L because of the potentially different activity 
factors for the participating BSs and the rest of the network. 
The number of active participating BSs interfering with the 
L lh BS is 12 = a i' which is a binomial random variable 

with probability mass function 

w £ {0, ... ,L — 1}. 

(3) 


/fi(w) = 


L- 1 

OJ 


D. Base station range distributions 

From Q, it is clear that the SINRs are dependent on the 
distances of the BSs from the origin, rather than the locations 
themselves. Thus, it is worthwhile to characterize these. Let 
Rk — ll^fcll and R-rn = ||a; m ||, where x m is now the m th 
farthest active BS. This is illustrated in Figure [l]for the closest 
active (Ri) and L th closest overall {Rl) BSs. The distribution 
of Rl is (36): 


Ir l 0 ) 


= -Attt - 2 2 ( A7rr2 ) £ 

r(L — 1)! ’ 


(4) 


Now, conditioned on 12 and the distance of the L th BS from 
origin, we present a useful lemma for understanding the 
distribution of the locations of the 12 active participating BSs. 

Lemma 1. Conditioned on the distance of the L ,h BS from the 
origin, Rl, the Cl active BSs closer to the origin than the L ,h 
BS are independent and identically distributed over b(o, Rl) 
with the location of each BS sampled uniformly at random 
from b(o, i?i). 

Proof: See Appendix [A] ■ 

Note that a point process is which a fixed number of points 
are distributed uniformly at random independent of each other 
in a given compact set is termed as a uniform Binomial Point 
Process (BPP) G3- Therefore, the above Lemma simply states 
that conditioned on Rl, Cl active BSs closer to the origin than 
the L th BS form a BPP on b(o, R L ). 

Using Lemma |T] we obtain the following distribution for the 
most dominant (i.e., closest active) interfere^ given the distance 
of the L th BS and 12. This is formally presented below. 

Lemma 2. The cumulative distribution function of the closest 
active BS distance ll\ given Rl and 12 is 

F RdR L ^ r \ R ^) = l-{^f^) ’ 0 < r < Rl- (5) 

Proof: See Appendix [B] ■ 

Let A = b(o, f?i)\b(o, f?i), where b (6,r) represents a ball 
of radius r centered at 6. We have the following lemma 
characterizing the distribution of the 12—1 remaining active 
BSs in the annular region A. 

Lemma 3. Conditioned on R\ and Rl, the 12—1 active 
BSs located inside the annular region b(o, i?i)\b(o, f?i) are 
distributed according to a uniform BPP. 

Proof: See Appendix [A] ■ 

III. Localization Performance 

As mentioned in Section [I] localization performance funda¬ 
mentally depends on the number of positioning observations, 
their accuracies, and the locations of the participating BSs 
relative to the device being localized. These three factors are 
all highly interdependent. For instance, both the number of BSs 
whose signals arrive with sufficiently-high SINRs to participate 
in the localization procedure and the quality of their positioning 
observations are strongly dependent on the interference field, 
which is itself driven by the realized network deployment. 
Thus, a full characterization of localization performance would 
take into account all possible geometric conditionings of the 
BSs over all possible channel conditions, a task which is 
extremely difficult. By taking into account even just the number 
of participating BSs averaged over the channel conditions, 
however, valuable insights into localization performance may 
be gleaned. 

A. Base station participation and geolocation performance 

A clear relationship exists between the number of BSs 
involved in positioning and a cellular network operator’s ability 
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Fig. 2. Impact of BS participation on meeting FCC E911 require¬ 
ments: Positioning performance of an OTDOA-like system with settings 
derived from |38) (see text). A performance curve must pass through both 
horizontal dashed lines in order to meet the FCC mandate. 


to meet its geolocation performance requirements (21) . In order 
to see this, consider handset-based TDOA positioning, such as 
Observed Time-Difference-Of-Arrival or OTDOA in LTE. The 
FCC E911 mandate requires that operators employing handset- 
based localization estimate the locations of mobile devices to 
within an accuracy of 50 m 67% of the time and 150 m 90% 
of the time (37). Using the model described in the previous 
section with a BS deployment density equivalent to that of 
an infinite hexagonal grid with 500 m intersite distances, a 
shadowing standard deviation of 8 dB, and a pathloss exponent 
of a = 3.76, an example localization system was simulated 
requiring a pre-processing SINR /3/7 = —13 dB for successful 
signal detection (38) . A fully-loaded network was assumed with 
no BS coordination (i.e., p = q = 1) and clock inaccuracies at 
each BS modeled normally with a standard deviation of 100 ns. 
The SINR threshold determines which BSs may successfully 
participate in the localization procedure, while an individual 
ranging observation is modeled as the sum of an exponentially- 
distributed non-line-of-sight bias with a mean of 30 m and a 
zero-mean Gaussian random variable with a variance based 
on the exact SINR and calculated using the well-known TOA 
ranging CRLB (39) using a signal bandwidth of 10 MHz. Lastly, 
positions were estimated using the algorithm presented in (5) 
with the strongest BS selected as the reference BS. The resulting 
localization performance is shown in Figure [2] for several 
minimum BS hearabilities. From this particular example, a 
network operator might conclude that, with universal frequency 
reuse on the positioning signals, a hearability of 6 BSs would 
be required to be in compliance with the requirements of the 
FCC E911 mandate. 

Motivated by this connection, we now formally propose 
a metric to study the number of BSs able to successfully 
participate in a localization procedure for a given target SINR 
/? and processing gain 7. 

Definition 1 (Participation metric). For a given BS deployment 
0 7 T, let T represent the maximum number of selectable BSs 


such that all successfully participate in a localization procedure. 
This can be mathematically defined as: 

1 / 

T = arg max l x 0 1 ( SINRfc(£) > 

1 k= 1 ^ 

Note that this metric differs from a traditional metric such as 
the CRLB in that it is not directly tied to a specific positioning 
accuracy value, though it is still closely related. Its advantage 
lies in its tractability when characterized over all the possible 
BS topologies, whereas the CRLB, on the other hand, does 
provide direct insight into achievable positioning accuracy, but 
only for deterministic networks, quickly becoming intractable 
to characterize over all BS topologies. 



B. Participation from a desired number of base stations 

While generally better localization performance can be 
attained by increasing BS participation, the probability of 
obtaining some desired number of successful BS connections is 
known to decrease sharply as the number increases. How this 
fundamental factor influencing localization system performance 
changes across different parameter sets is exactly the objective 
of our study. Thus, it is desirable to understand exactly how T 
is impacted by the network design parameters and propagation 
effects, averaged over the geometric conditions and channel 
variations. 


Definition 2 (L-localizability probability). For a given 0 7 <1>, 
a mobile device is said to be L-localizable if at least L BSs 
may successfully participate in the localization procedure. The 
probability of this occurring is simply: 


Pl = P (T > L) = E 


L 


n* 


^SINR fc (L) > 



(7) 


It may be insightful to step back and consider a specific 
application of P L . Note that the first objective in any location 
system is to make sure that the device to be located can detect 
positioning signals from a sufficient number of BSs (i.e., that 
it is localizable). By this, we simply mean that an estimate of 
the device’s location can be found without ambiguity. In the 
noiseless case, this means that there can only be one solution. In 
the noisy case, this means that there is a single global minimum 
to the appropriate cost function. Commonly-accepted minimum 
values of L for the unambiguous operation of a localization 
system in the R 2 plane are 2, 3, and 4 for AOA, TOA/RSS, and 
TDOA, respectively. Thus, for example, P 4 can be equivalently 
thought of as the coverage probability of TDOA positioning. 

For two edge cases of interest, that with no BS coordination 
(p = 1) and that with perfect BS coordination (p = 0), it is 
straightforward to infer from 0 that 


1 (SINR fe (L) >/?)>! (SINR;(L) > j3) (8) 


for all h < l < L. This simply means that the received SINR 
from a BS farther from the mobile device is lower than that of 
a closer BS, implying that the probability of A-localizabi Iity 
in 0 can be equivalently expressed as 

P L = E[l(SINR t (L) >0/7)]. (9) 
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Fig. 3. Lower BOUNDING THE INTERFERENCE: As in Figure JT] trans¬ 
mitting/quiet BSs are denoted by filled/hollow circles and the device to be 
localized is denoted by the filled square. 


With partial BS coordination (0 < p < 1), the relationship in 
([ 8 ]) does not hold in certain corner cases. These rare occurrences 
have a minimal impact on our analysis, and we proceed by using 
the expression for P L in (|9]> for all p. In the numerical results, 
we will validate that 0 is in fact an accurate approximation 
of P L for all p by comparing with the true P L , which jointly 
considers the SINRs of all participating BSs. 


C. A simple bound on P L 


all participating BSs transmit simultaneously, the above bound 
reduces to a simple closed-form expression. 

Corollary 1.1 (The special case of p = 1). When all 
participating BSs transmit simultaneously, 

C _ 2 \ L —1 

!- (j- (L - 2 )) J > (U) 

for p < jyfx, an d zero otherwise. 

Proof: When p = 1, fn(L — 1) = 1 and the result follows 
from JTO) . ■ 

From this result, we see that as greater numbers of participating 
BSs are desired, the probability of attaining those numbers of 
BSs with sufficiently strong connections decreases dramatically. 
Note that by rearranging CD- we also obtain a useful lower 
bound on the processing gain required to achieve some desired 
L-localizabity probability. 

Corollary 1.2 (Lower bound on processing gain with no BS 
coordination). The processing gain required to reach a desired 
L-coverage probability P L when p = 1 is lower-bounded by 

7 >^(l-P L ^)” 7 +L-2^j . (12) 


D. Approximations of P L 
Let 

n 

i 1 = Y J pm\~ a (i3) 

i=2 


We now derive an upper bound on L-coverage probability 
by (i) considering the closest interferer’s location exactly, (ii) 
placing the remaining 12—1 active interior interferes at the edge 
of the field of participating BSs, and (iii) ignoring thermal noise 
and all distant BSs. We illustrate this setup in Figure [3] which 
applies the above considerations to the network realization of 
Figure [Tj Note that the interference from all BSs beyond the 
L th BS is ignored, while the active interior BSs have been 
pushed out to the same distance as the L th BS. The dominant 
interferer, however, has been left untouched. 

Theorem 1 (Upper bound on L-coverage probability). The 
probability that a device can successfully utilize a desired L 
BSs for positioning is upper bounded by 

P L (p,g,a,/3,7, A) < ^ ^1 - ^ - (w - l)j j /nM, 

( 10 ) 

where \ = min{L — 1, fy/P J}. 

Proof: See Appendix [C] ■ 

One thing we notice immediately from ( fTO) is that the bound 
is independent from the BS deployment density A. This is 
in line with similar observations made in coverage analyses 
for interference-limited networks (i.e., network interference 
dominates thermal noise) (25), (26) . We will return to this 
observation at the end of this section. In the special case when 


be the aggregate interference due to the active BSs (if any) 
between X\ and x /. Furthermore, let 

OO 

* 2 = £ b 5 P\\ Xj \\- a (14) 

j=L + l 

be the aggregate interference due to the infinite field of BSs 
located further than the L th BS. Clearly then, when O > 1, 


SINR l (L) 


PRf a 


PR L a 

+ li + I 2 + c 2 


(15) 


In this section, we will approximate (15) and in turn P L by 
making the following assumption. 


Assumption 1. We assume that if the dominant interferer 
(i.e., the closest active BS) is considered exactly, the remaining 
interference terms in CD- namely X\ and X 2 , may be accurately 
approximated by their means conditioned on Ri, R^, and O. 


This type of dominant interferer analysis has been employed 
previously with desirable results (27), (28) and will yield 
remarkably simple, yet accurate, approximations in our analysis, 
as will be demonstrated later. This is because considering If 
exactly and 1\ using its mean (conditioned on If ), results in 
an accurate approximation of the total interference due to the 
BSs closer than xl, which is typically the performance-limiting 
term. 
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In interference-limited networks. Assumption |T] allows us to 
replace SINRi(L) in (|9]» by 


SIRl(L) 


PR~ L a 


PRf a + E[l 1 \R ll R L ,n] + E[l 2 \R L \' 


(16) 


where SIR stands for the signal-to-interference ratio. As the 
deployment density grows, the assumption of interference- 
limited networks increases in validity, and since our work 
is motivated by the reuse of existing infrastructures, such 
as cellular networks with small cell extensions which have 
relatively high deployment densities, we aver that it is quite 
reasonable. Expressions for the above conditional means are 
presented in the following lemmas. 


Lemma 4. The expected value of T\ conditioned on R \, Rl, 
and Q is 


E[h\R x ,R L ,Sl) 


2 P(fl - 1) R 2 L - a - R\- a 

2 — 0! ' R 2 l -R\ 


(17) 


for fl > 1 , and zero otherwise. 

Proof: See Appendix [D] ■ 

Lemma 5. The expected value of !■> conditioned on R^ is 


E[I 2 | Rl\ 


2PnqX R 2 ~ a 
o-2 L 


(18) 


for a > 2, and unbounded otherwise. 

Proof: See Appendix |E] ■ 

By inserting the results of Lemmas [4] and [5] into ( p~6l >, we obtain 
increasingly specialized (and simpler) approximations of the 
L-localizability probability. In order to proceed, however, we 
must first account for the case when il = 0, in which case R x in 
© has no meaning. Thus, we first derive an approximation of 
P L for the special 11 = 0 case, which incidentally corresponds 
to the perfect coordination (p = 0 ) scenario. 


Proposition 1. Under Assumption [7] the probability of L- 
localizability with perfect BS coordination in interference- 
limited netw’orks is 

L—t 2 ( tt ~ 2 V 

P L (0,g,a,/3,7, A) = 1 - ' , (19) 

e=o 

Proof: See Appendix [F] ■ 

As noted in the proof of above Proposition, (19) can be 
interpreted as the probability that at least L BSs (modeled by 
a PPP) lie inside b ^o, ^ j • Moreover, note that the 

density of the BS deployment does not appear to affect P L . 
With this special case out of the way, we now proceed, starting 
with our most general result. 


Theorem 2. Under Assumption [7J the probability that a mobile 
device is able to localize itself using L BSs in interference- 
limited networks is P L (p, q , ct, p, 7 , A) = 


L—l 


/ 

1 ^ ^ £ 2<z/3/T' 

^ £=0 


(S=2.Y 

J\ 


/n( 0 ) 


4(A7t) l 


L—l 




/»oo 

/ 1 

r J 

'0 J 

» \ 


2(ui-l) 
2-a 


. ~ r i 

,,2 — f 2 

' r 1 l 


J 2nq\ 2-a 


y-2 ' L 


h(' 


,2 _ ~2\ui-l 2(L-ui)-l \ vr 


' 1) 


uie~ L dfi dr L . 



( 20 ) 


Proof: See Appendix |G] ■ 

Thus, we obtain an expression which requires only double 
integration in order to determine P L , and though the outside 
integral has an infinite integration bound, the expression is 
not difficult to evaluate numerically thanks to the decaying 
exponential term in the integrand. To appreciate the value 
of this approximation, consider the p = 1 scenario in which 
case the first term in ( [ 20 ] ) disappears and the summand need 
only be evaluated at w = L — 1. Though the authors did 
not have localization in mind, the exact expression for the 
fc-coverage problem in the absence of fast fading in |40] Corr. 
7] applies directly to P L for p = q = 1. A cursory glance at 
the exact expression in [401 reveals its cumbersome nature, 
involving sums over many multi-fold integrals even in the 
absence of thermal noise, which is considerably more complex 
in comparison to ( |20| . Moreover, by applying a slight secondary 
approximation and multiplying E[I 2 | 7 ?l] in ( fT 6 | > by E [Rf\/R\, 
where E[i?f] = [ 361, ( | 20 | > is further reduced, yielding a 

single-integral expression applicable to the general case. 


Corollary 2.1. The general expression for P L in ( |20| ) may 
be further (and consequently, less reliably) approximated for 
interference-limited networks by: P L (p, q, a, ft, 7 , A) = 


L—l 

Pl(0, q , a, P, 7 , A)/q( 0) + 2^ /n(w)w 

0J= 1 


1 



2 (w — 1) x 2 — x a 
2 — a x 2 — 1 


2ga! 2 < 7\ 

a —2 " P) 



0J — 1 

-dx. ( 21 ) 


Proof: See Appendix [H] ■ 

This expression is instantly and accurately evaluated in a 
numerical computation environment such as MATLAB. We now 
briefly describe how the introduction of the additional terms 
into ( fl 6 | i reduces Q. The key is that after some algebraic 
manipulation, © may now be expressed as follows (see © 
in the appendix): SIRl(L) ~ 



1 


2(0-1) ( Hi ) ( Hi ) 1 2q ( R l \~ 

2 —q - (Rl) 2 _ 1 ^ a-2 \r 1 J 


( 22 ) 


which contains i?i and llr only through their ratio X = 
Rl/R\. The cumulative distribution function of this ratio 
admits the simple form of Fx{x) = (1 — l/a; 2 )^ 2 , which 
is readily obtained from ([5]). For the special case of a - 4, 
this additional approximation needn’t be used, however, as a 
single integral expression is obtained directly from ( | 20 | >. 

Corollary 2.2. Under Assumption [7] and for the special case 
of a = 4 in interference-limited networks, 

x 

Pl(p, 9,4, p, 7 , A) = ^2 /n(w) 

cd —0 






























1 


f RL (r)dr, (23) 


IV. Numerical Results and Discussion 



\Jl/P- TtqXr 2 + ^ 4 ^ 


where \ = min{L — 1 , Lt//3J }• 

Proof: See Appendix [1] ■ 

For fully-loaded networks (i.e., p = q = 1), ( |23[ ) simplifies 
further to the following expression. 

Corollary 2.3. Under Assumption [7] and for the special case 
of a = 4 in fully-loaded interference-limited networks, 


P L (1> 1,4, /3,7, A) 




-V/P-(L- 1) 

7T X 



1 

yJl/P- nXr 2 + (L 4 2)2 



L —1 

(24) 


/or /3 < jffr and zero otherwise. 


For larger values of L and higher values of p, it may be 
reasonable to ignore the I 2 term in © under the assumption 
that the interference due to the f l active participating BSs 
dominates that of the rest of the network. By doing so, 
expressions not requiring any integration may be found for the 
special cases covered in Corollaries |2.2| and |2.3| 

Corollary 2.4. Under Assumption [7] a = 4, P L for networks 
limited by the interference from the nearby BSs participating 
in the localization procedure is P L (p, q, 4, /3, 7 , A) = 


Y ( 1 


uj—0 



where \ = min{L — 1, [_7//3J }• 


(25) 


Proof: Consider the limit of ( |23j ) as q —> 0. The 
parenthetical term becomes independent of r and may be 
pulled out of the integral, while the remaining integral covers 
the probability density function fn L (r) over the entirety of its 
domain, thus evaluating to 1 . ■ 


Corollary 2.5. Under Assumption [7] and when a- = 4, P L 
for fully-loaded netw’orks limited by the interference from the 
nearby BSs participating in the localization procedure is 


P L (1,1,4,/3, 7 ,A) 



1 

\J l/P + 


for /3 < jyzj and zero otherwise. 



(26) 


Interestingly, when interference from the participating BSs is 
ignored in ( fl9| and when distant network interference is ignored 
in ( |25j ) and ( |26[ ), the resulting expressions are independent of 
the BS deployment density. This is termed as scale invariance 
in the literature and is known to hold for most of the metrics 
that are derived from SINR when the thermal noise is negligible 
or ignored (interference-limited scenarios), e.g., see |[25), [ 26) . 


We now study the tightness of the bound and the approxi¬ 
mations, as well as gather insights from numerical results. For 
consistency, all results were gathered using a BS deployment 
density equivalent to that of an infinite hexagonal grid with 
500 m intersite distances (i.e., A = 2/(v / 3 • 500 2 m 2 )) and a 
shadowing standard deviation of 8 dB. While our analytical 
expressions were derived by focusing solely on the post¬ 
processing SINR from the L th BS surpassing /3, the truth 
data, which is compiled through simulation using an average 
of 1000 BSs, contains no such approximation and considers the 
SINRs from all L participating BSs jointly with each required 
to surpass /?. Let us first consider P L for L = 4 in the case 
of a = 4 and fully-loaded networks (p = q = 1), the setup 
considered in Figure [4] To interpret the meaning of this figure, 
it may be helpful to provide a simple example. Consider a 
cellular localization system employing TDOA observations 
and requiring a post-processing SINR of at least —6 dB for 
successful operation after experiencing a 10 dB processing 
gain due to its integration time (i.e. /J /7 = —16 dB). Figure [f] 
tells us that when all BSs in the network are fully-loaded (i.e., 
actively transmitting), unambiguous TDOA positioning (i.e., 
using L = 4 or more BSs) is possible only 50% of the time, a 
result consistent with the fact that universal frequency reuse 
(e.g., using synchronization signals) was deemed untenable for 
positioning signals in LTE. Also in this figure, the truth data 
is compared against the closed-form bound in Corollary 0 
the nice closed-form approximation in Corollary |2.5| which 
ignores distant interference, as well as the simple single-integral 
approximation in Corollary |2.3| We note that the bound is tight 
for the high-reliability regime, which is the range of greatest 
interest (localizability probabilities greater than 0.75). For low 
probabilities of localizability, the bound diverges by a maximum 
of around 3 dB. However, this low coverage range is not very 
desirable for system design, anyway. In addition to the bound, 
the closed-form approximation performs well for P L > 0.6, 
while the single-integral approximation provides a remarkably 
tight approximation which is invariably indistinguishable from 
truth across all SINR values. The key factor in enabling these 
accurate approximations is that the locations of the strongest 
BS (i.e., the dominant interferer) and the L th BS are considered 
exactly. 

Inspired by the results in Figure |4l we now delve deeper 
into the approximation of Theorem [2] in ( |20) and its single¬ 
integral counterpart in ( | 2 Tj ) to examine how they perform 
under more general conditions, such as across the range of 
realistic pathloss exponents. Figure [5] shows that the former 
remains almost indistinguishable from truth for all a values, 
while the latter is very tight in the high-reliability regime 
for a > 3 and exhibits a noticeable deviation (amounting 
to less than 1 dB) when a = 3. For P L < 0.5, the single¬ 
integral approximation deviates significantly, however, due 
to the fact that the interference from beyond the if' BS 
becomes more influential with a higher SINR requirement. This 
deviation is further accentuated by our intentional selection 
of disadvantageous p and q values, where p = 2/2, partially 
de-emphasizes the interference from the closest BSs and < 7=1 
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Fig. 4. An INITIAL COMPARISON: This figure compares the true Pl with the 
bound and several approximations for L = 4, a = 4, and p = q = 1. The 
most accurate approximation is indistinguishable from truth, and requires a 
simple analytical evaluation instead of a tedious Monte Carlo approach. 


is the worst case q for the single-integral approximation since it 
maximally emphasizes the distant interference term where the 
additional approximation lies. Nevertheless, this single-integral 
approximation for the general parameter case performs well 
even under these adverse conditions for the higher P L values 
of interest, and will only improve with more favorable values 
of p and q. On the other hand, subsequent figures will make 
it evident that the approximation in Theorem [2] and its exact 
derivatives remain extremely accurate across the ranges of all 
network parameter values. 


As discussed in Corollary 1.2 the upper-bound on coverage 


provides a useful closed-form lower-bound on the processing 
gain required to obtain a desired A-localizabi Iity probability, 
which is likely a specification set prior to system design. From 
Figure[6] it is clear that for p = 1 and high desired values of P L , 
the bound provides a tight criterion for the minimum processing 
gain required. Unfortunately, a closer look reveals that the 
processing gains required are unpractically high, highlighting 
the need for some sort of interference mitigation. Moreover, 
consider that both bounds are very tight in the region of 
interest despite the fact that the interference from beyond 
the L th BS is fully present (i.e., q = 1), while the bound 
derivation essentially assumes q = 0. This provides more 
insight into the nature of the interference, revealing that the 
localization performance is limited by the interference from 
other participating BSs (which is consistent with conclusions 
drawn from LTE positioning simulations |21|), prompting 
one to consider coordination amongst the BSs as the desired 
interference mitigation technique. This becomes even more 
evident as the propagation conditions become more severe, 
such as indoors, where higher path loss exponents demand 
higher processing gains to achieve the same P L as seen in 
Figure [7] 

Through Figure [8] we demonstrate the value of mitigating 
interference through BS coordination. With the BSs beyond 
the L th BS fully-loaded, we note an increasing benefit from BS 



Pre-processing SINR, p/y (dB) 
(a) Theorem [5] 



Fig. 5. General Theorem|2]approximations across a range of pathloss 
exponents for L = 4, p = 2/3, and q = 1. 



-25 -20 -15 -10 

Post-processing SINR required, p (dB) 


Fig. 6. The importance of interference mitigation: Processing gains 
required to achieve Pl = 0.8 without any sort of interference mitigation 
(p = 1) for q = 1 and a = 4. 
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Fig. 7. THE IMPACT OF path LOSS: Processing gains required for desired 
Pl values grow exponentially with the pathloss exponent. (Pl = 0.8, L = 4, 

p = q = i) 



Fig. 9. The cost OF INVOLVING more BSs: This figure underscores the 
high cost of obtaining a greater desired number of BSs, particularly in the 
high-reliability regime, (a = 4, p = 1/2, q = 3/4.) 



Fig. 8. THE BENEFIT OF BS COORDINATION: These L-localizability 
probabilities emphasize the gains achievable through the coordinated mitigation 
of participating BS interference. (L = 4, a = 4, q = 1.) 


coordination as the degree of the coordination is incrementally 
increased. The large improvement in hearability from p = 1 to 
p = 0 again sheds light on the issue of interference from nearby 
BSs. In fact, at its least detrimental point, the interference from 
the L — 1 closest BSs accounts for nearly a 10 dB drop in 
signal quality. That’s roughly the equivalent of going from 
successfully obtaining eleven to obtaining only four BSs in the 
network considered in Figure [9] While coordination is able to 
help significantly, a closer look at Figures [8] and [9] makes it 
clear that more is necessary in order to reach the truly high 
L-localizability probabilities at reasonable SINR values. The 
most prevalent technique to further reduce the interference is 
that of thinning the interference field by separating nearby 
transmitters in frequency, which is discussed next. 


A. Application to frequency reuse 

An effective technique to reduce interference from co-located 
transmitters is that of frequency reuse, commonly included in 
wireless standards [41 ]. If a total of K frequency bands are 
available and we independently assign one of the bands to 
each a; £ $ with equal probability, we can easily incorporate 
frequency reuse into our model by considering the transmission 
activity on each band separately using independent PPPs whose 
densities are that of the original PPP thinned by the frequency 
reuse factor K. Let us denote the number of detectable BSs 
in the fc th band by 77. For L-localizability, we thus require 
ICitLi n k L L. The main technical difference between this 
analysis and that already presented is that we are now obliged 
to evaluate the probability of being able to detect exactly n 
BSs in a given frequency band. Under the localization setup we 
consider, it is not difficult to show that when p = q this can be 
readily determined by considering the difference between the 
probabilities of n and n +1 localizability. The p = q condition 
isn’t particularly restrictive, as it still allows us to model average 
network load and only removes the ability to perform per- 
band BS coordination (which isn’t done in practice, anyway). 
From here, the expression for the L-localizability probability 
with frequency reuse is procured by simply considering all 
combinations of {ni ,..., uk} for which Y^k=i n fc > L. This 
is formally stated in the following theorem. 


Theorem 3 (L-localizability probability with random frequency 
reuse). The L-localizability probability with random frequency 
reuse and reuse factor K is: 


L—l 


K 


PE(p,«,a,/3,7. A) = 1 - ^2 FI 


/—0 


K 


P«i P,<?,a,/3,7, — -P„ i+ i [p, q, a, /3,7, — 


K 


(27) 


This expression is easily implemented in software using a 
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Fig. 10. THE value OF FREQUENCY REUSE is undeniable for improving P L . 
Here, L = 4, a = 4, p = q = 1. 



Fig. 11. CONVERGENCE TO PPP: Hearability results gathered using a regular 
grid model converge to those gathered using a PPP model as the shadowing 
standard deviation a a increases. (3f'y = —10 dB, a = 4, p = q = 1.) 


recursive function and any of the previous expressions for 
evaluating P L may be substituted into ( |27[ i. Using the worst- 
case interference scenario (p = q = 1), the value of frequency 
reuse is clearly demonstrated in Figure [TO] for K = 3 and 6, 
which correspond to the values used in LTE when positioning 
using cell-specific reference signals (CRS) and positioning 
reference signals (PRS), respectively. Since random frequency 
reuse makes no attempt to intelligently separate interference 
sources, gains from planned frequency reuse are likely to be 
even better. 


B. Example LTE use case 

Although the model presented in this paper is general and 
not tied to any specific localization system, here we provide an 
example of how our model may be initially related to an existing 
technique. Consider the scenario of OTDOA positioning in 
LTE using the PRS, for which it is reasonable to assume that 
up to six close BSs may be selected to localize a mobile 
device without interfering with each other due to a frequency 
reuse factor of six among the PRS. BSs farther away may still 
interfere, however. This setup may then be treated as one with 
perfect coordination (p = 0) among the L = 6 participating 
BSs, while the expected interference from beyond the sixth 
BS is adjusted using q < 1 in accordance with the average 
network load. The processing gain 7 ss 15 dB from the length- 
31 Gold sequences used in the PRS while a = 3.76 per 3GPP 


assumptions [20]. Furthermore, even PRS power boosting may 
be accounted for in a preliminary way by considering the power 
boost at the participating BSs as an equivalent power reduction 
at the interfering BSs. Practically, this may be implemented by 
reducing the value of q, the ultimate value of which will depend 
on the amount the participating BSs’ powers are boosted. 


C. Similarity between PPP and hexagonal grid hearability 

We now briefly revisit our assumption of BSs distributed 
according to a PPP and show that this model yields hearability 



8 9 10 11 12 13 14 15 16 17 18 19 20 

Number of BSs, L 


Fig. 12. Comparison with Frequency Reuse: When K > 1, the models 
agree for greater values of L even at a lower cr s = 8 (111. A required L > 10 
is not typically expected, i/3/7 = —10 dB, a = 4, p = q = 1, K = 6.) 


results which are comparable to those yielded by the well- 
accepted hexagonal grid model. As discussed in |[32)-[34), 
increased shadowing causes more regular network models to 
appear Poisson when considering received signal strength-based 
metrics. In Figure [TT] hearability results using the two models 
are compared for several increasing shadowing strengths (note 
that, as discussed previously, the PPP results are invariant to 
shadowing strength), p = q = 1, and /d/'y = —10 dB, and it is 
clear that the models are converging. Furthermore, frequency 
reuse causes the models to provide similar results for higher 
values of L even at lower shadowing strengths. This is shown 
in Figure [12] for K = G, <j s = 8 dB, and the other parameters 
the same as above. While there will always be some difference 
in the hearability results using the two models, the results using 
actual deployments likely lie somewhere in between the two, 
as was discovered in |[25j for cellular coverage. 
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V. Conclusion 

We have employed concepts from point process theory and 
stochastic geometry in order to provide a new tractable model 
for studying localization in cellular networks. The model is 
accompanied by an analysis of hearability, which is shown 
to be an important metric providing insight into fundamental 
localization performance, and easy-to-use analytical expressions 
for this metric are provided. This is in contrast to most 
previous approaches, which either provide insights specific 
to deterministic deployments or rely on time-consuming 
simulations, neither of which allow for the drawing of general 
insights. While more base stations participating in localization 
is beneficial to its accuracy, the analytical results show that 
obtaining an increasing number of successful base stations 
connections decreases exponentially. The primary culprit for 
this is the interference due to the participating base stations 
themselves. When base stations coordinate (e.g., through a 
network controller), more of them are able to successfully 
participate in the localization procedure. However, in the end, 
it is clear that the greatest gains in hearability are realized 
through frequency reuse. This is consistent with the findings 
during the standardization of LTE positioning, where this 
conclusion was reached after performing many complex system- 
level simulations. 

We have only scratched the surface of the model, which 
may be used to extend this work in many different ways. 
Besides hearability, understanding how other metrics, such as 
ones focused on network geometry, are affected by changes 
in design parameters would yield valuable contributions to the 
localization literature. In addition, recognizing how directional 
antennas and power control affect localization performance 
would be helpful to cellular system designers. Lastly, location 
information can be gathered in diverse ways, for example by 
simply knowing that one is able to hear a low-power femtocell, 
which lends itself nicely to the consideration of different 
classes of base stations through multi-tiered heterogeneous 
networks [26] . 


Appendix 

A. Proof of Lemmas |7] andjJj 

The proof strategy used here is quite standard, e.g., see 1 12 
Theorem 2.9]. For the proof, let i < j, n < N, K be 
the number of BSs between Xi and Xj, A = b(o, ||a:j|| + 
<5/2)\b(o, \\xi\\-8/2), B = b(o, ||a: i ||+e/2)\b(o, \\xj\\-e/2), 
C = b(o, ll^ll — e/2)\b(o, Hscill + 5/2), D C C, b(0,r) 
represent a ball of radius r centered at 6 , 5 and e be small 
modifications to the regions of A, B, and C, and Ng be the 
number of points in region PL. Then, P(Ng = n\xi,xf) = 


(\\D\) n e~ x \ D \/n\ (A(|C| - |L>|))*- n e - A (l c 'H^I) 


(A|C|) N e -A l c l /K! 


(K — n)\ 


\D[ 

\C\ 


1 - 


\D _| 
\C\ 


H-r. 


which shows that the H BSs inside the annular region 
lim^e^o C = b(o, ||Xj||)\b(o, ||xi||) are distributed uniformly 
at random independent of each other. In other words, they form 
a uniform Binomial Point Process (BPP) over the annular region 
b(o, ||xj||)\b(o, Hsill). By letting i = 0, j = L, and defining 
Xq — o (i.e., ||xo|| = 0), the BSs inside the circular region 
b(o, ||xl||) are shown to form a uniform BPP, thus proving 
Lemma [T] By letting i = 1 and j = L, Lemma [3] is proved. 


B. Proof of Lemma [2] 

Let x be an active BS that is uniformly distributed inside 
the circular region b(o, R /,). Then, for il > 1, 

= 1 - P(Ri > r\R L ,n) 

= 1 - P(min{x : ||x|| < R L } > r\R L} Ll) 

( = } 1- [] P (ll*ll > r \ R L) 

{x-.\\x\\<R l } 


(tO 


=’i- n 


7r R 2 r — nr 2 


7 tR 2 , 

{x-.\\x\\<R l } 


= 1 - 


Rj-r* 

R 2 l 


where (a) and (6) follow from Lemma [I] 

C. Proof of Theorem [7] 

For fl > 1, 


P L = P 


< P 


PR~L a 


> - 




^ p 


(b) 

< P 


,Ei=! PR- a -PR- L 
PRl a 

X?=1 PR- a -PR- L a “ 

f _S 

[Rf a + zL K a 'y J 

( Rf a ^ 

R^ + (n-i)R- L a - 7 y 


= P At > Rl (i - {n - i) 


- p ( min{||x||: \\x\\ < R L } > Rl ( ^ ~ H - 1) 


0 |p 

= bn 




1 ( - < H _1 

,7 


lim P(N d = n\N c = N B = 1, N A = 1) 

<5,e->0 

= P (N D = n, N c =*,N b = 1, N a = 1) 

S,e —>0 P (Nc = Ng = 1, N A = 1) 

P(jVp = n, N C \p = X-n,N B = l,N A = l) 
5,e—y 0 P(iV C = j = 1, N A = 1) 

_ P(N d = n) P(N C \d = H - n) 

” P(^c = «) 


where (a) follows from Rn+i = Rl and pulling the dominant 
interferer out of the sum, (b) from the fact that R, < Rl for 
i < H, (c) since R\ is the smallest among all interferers, and 
(d) from the fact that given Rl, all the BSs are uniformly dis¬ 
tributed in the circle of radius Rl, followed by some algebraic 
manipulations. The final result is obtained by deconditioning 
over Q > 1, resulting in an expression which also trivially 
valid for the fl = 0 case. 
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D. Proof of Lemma |4] 

Let x be an active BS inside the annular region 
b(o, i?i)\[b(o, i?i). Then, 

l-Rl 


E[P\\x\\- a \Ri,R L } = I Pr~ 

J Ri R l - R{ 

rR-L 

r 1_a dr= 


2 r 


dr 


2 P 


2 P 


Rl-R\ 


R\ -R\ \ 2 -a 


Rl 


r=Ri , 


2 P R,f n - /(, 


<2—o: 


2-a 


jd 2 r>2 

U L R i 


where (a) follows from the fact that x is uniformly dis¬ 
tributed inside the annular region b(o, f?i)\b(o, Ri) (proved 
in Lemma |3j. Lastly, the mean of the sum of the interfer¬ 
ence from all interferes located inside the annular region 
b(o, f?i)\b(o, Ri) is simply the sum of their individual means, 
from which the result follows. 


E. Proof of Lemma [5] 

In order to accommodate partial network loading, q is the 
fraction of the far-off BSs transmitting throughout the local¬ 
ization procedure, which can be modeled using an equivalent 
thinned PPP ip with density qX. For a > 2, EfZ^I Rl\ = 


E p IMI 

p2n poo 


^qX 


PM 


*dx 


lb <={o,R L ) 


= PqX 


r a r dr dd = 2PnqX 


r 1_a dr, 


>Rl 


’Rl 


where (a) follows from Campbell’s theorem Jl2) , and the result 
follows by solving the final integral (assuming a > 2). 


F. Proof of Proposition [7] 

The L-localizability probability in this case is P L |n =0 = 


-Rl 


= E 


( PR 


— a. 




Rl 


VE[X 2 |7? l ] ~ 1 )\ 

( a-2 


= E 


Rl 


1 , T > 

- o-2 L 


\2TrqX(3/'y 


>R 2 l 
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L—l 


! a—2 \t 


(a) , v—r_ 

= 1 — \ e 2< ?/ 3 /7 ■ 




t\ 


where (a) follows from the fact that the previous expression is 
simply the probability that at least L BSs lie inside the region 


b (°V 


*~2 ) 
2ttijA/3/7 ) ■ 


G. Proof of Theorem [2] 

The general L-localizability probability is 


Pl = P ( SIR l > 


= Er 


-Rl 


-Ri 


1 ( SIR l > - 

7 


Ri,R L ,n 


R Ll V 


(28) 


The SIRl term in the above expression is 


SIRl = 


_ PJRT _ 

PRf a + Ep 1 \R 1 ,R L ,Xl] + E[X 2 | R l \ 

PR~ L a 


DD _ c* i 2P(Q—1) R l R ± , 2P-Kq\ d2 — ol 

rn l + 2—a • Rl-R\ + a—2 U L 


(29) 

(30) 


The result follows by substituting this expression for SIRl 
back in |28|, followed by integrating over the densities of R t , 
Rl, and fl > 1, and treating the fi = 0 case separately using 
Proposition [I] 


H. Proof of Corollary 2.1 


As in the Proof of Theorem [2] above, P L can be expressed 
in terms of SIRl as ( [28] ), where SIRl is given by ( |30| >. In 
this Corollary, we multiply the E[l21 Rl] term of SIRl by 
E [Rf\/R\ = 1/(7tA R\) to get a simpler approximation for 
SIRl, and consequently for P L , as follows: 


SIRl = 


PRY 


PRf a 


2P(fi-l) 
' 2—CL 


R,7 a -Rf 

d 2 p2 

K l PC 1 


2Pirq\ r/2-a E[P?] 
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2 q 
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(31) 


where (a) follows by substituting IiiJR\ — > X. The distri¬ 
bution of X can be readily obtained from |5]i as Fx (x) = 
(1 — l/x 2 ) n . Note that after this substitution, the dependence 
of SIRl on Ri and Rl is completely captured through X. The 
final expression is now obtained by substituting © in ( [28] ), 
followed by integrating over the support of X and subsequently 
summing over the probability mass of XI, excluding XI - 0 

which is treated separately using Proposition [I] The key here is 
that instead of the double integration over Ri and Rl, single 
integration over their ratio X suffices in this case. 


I. Proof of Corollary |2.2| 

As in the proof of Corollary |2.1[ we start with the expression 
of P L given by ( [28] ) in the proof of Theorem [2] The SIRl term 
appearing in the expression of P L is given by ( [30] ). In this 
Corollary, we simplify the SIRl expression for the case of 
a = 4. Starting with the SIRl expression given by ([30]): 


SIRl =- 

PRf a 


RRl a 

2P(f2-l) R 2 , q -P? “ 
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(32) 


where (a) follows by defining X = Sf- and substituting a = 4. 

^1 i . . 

Using ( |32p , we now simplify SIRl > /?/7 which will then be 
used to derive P L from ( [28] ). Defining Y = X 2 we have 

SIRl > (3/7 W V 2 + (XI - 1 )Y < kT 1 
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(n-i) 


/ -i , (^-i ) 2 
< K H--- 
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(c) 


0 < Y < \ k 1 + 


, (»~ 1) 2 (»~ 1 ) 


X 2 < a/k" 1 
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^ 1 < X < 




fp + ^ ^ 


u + 4 2 

< -Ri < iiz,, 


where n 1 = 7//3 — tt\R\ in (a), (6) follows from the fact 
that Y > 1, and (c) from Y = X 2 . Note that (6) and (c) 
require « _1 > f " j . Step (d) follows from X > 1. The 
earlier condition on re -1 is replaced by a more strict condition 
k~ 1 > ft in (d). Substituting this back in ( [28] ) and doing some 
algebraic manipulations, we get 


Pl = E n 


= E n 



f R\- 

Rl \ 

n- 

- 


7«- i + (n - i)2 ( v> 




R2 L 

\ J 





from which the result follows by simply deconditioning over ft. 
Due to k _ 1 > ft, the integration limits are from 0 to J . 
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